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Abstract 

We have studied quantum phase transition induced by a quench in different one dimensional spin 
systems. Our analysis is based on the dynamical mechanism which envisages nonadiabaticity in 
the vicinity of the critical point. This causes spin fluctuation which leads to the random fluctuation 
of the Berry phase factor acquired by a spin state when the ground state of the system evolves in 
a closed path. The two-point correlation of this phase factor is associated with the probability of 
the formation of defects. In this framework, we have estimated the density of defects produced in 
several one dimensional spin chains. At the critical region, the entanglement entropy of a block 
of L spins with the rest of the system is also estimated which is found to increase logarithmically 
with L. The dependence on the quench time puts a constraint on the block size L. It is also 
pointed out that the Lipkin-Meshkov- Click model in point-splitting regularized form appears as a 
combination of the XXX model and Ising model with magnetic field in the negative z-axis. This 
unveils the underlying conformal symmetry at criticality which is lost in the sharp point limit. 
Our analysis shows that the density of defects as well as the scaling behavior of the entanglement 
entropy follows a universal behavior in all these systems. 
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I. INTRODUCTION 

Kibble [l| has pointed out that in second order phase transitions if the system is driven at 
a fixed rate characterized by a quench time its evolution cannot be adiabatic close to the 
critical point. The nonadiabatic evolution in the critical region produces defects such that 
the system becomes a mosaic of ordered domains whose finite size depends on the transition 
rate. In the cosmological scenario Kibble considered relativistic causality to set the size of 
the domains. Later, Zurek ^ suggested the dynamical mechanism based on the universality 
of critical slowing down and predicted that the size of the ordered domains scales with the 
transition time Tg as where x is a critical exponent. The Kibble-Zurek mechanism (KZM) 
in second order phase transition essentially involves thermal fluctuation which initiates sym- 
metry breaking. In recent times, the zero temperature quantum phase transition (QPT) has 
been studied in the light of KZM and the formation of defects in one dimensional transverse 
Ising model has been investigated in details by several authors 3|-[5|. In QPT of the spin 
chains, apart from the study of the formation of defects it is also important to examine 
how entangled various parts of the system are with each other. Recently, it has been shown 
that for a quantum spin system in a mixed state the entanglement of formation of a pair 
of nearest neighbor spins i.e. concurrence is related to the Berry phase factor acquired by 
a spin state when the ground state of the system evolves in a closed path {g-S]. When 
a bipartite quantum system is in a pure state the measure of entanglement between two 
subsystems is given by Von-Neumann entropy. When the bipartite system is in a mixed 
state the entanglement of formation given byconcurrence has the property that it reduces 
to the Von-Neumann entropy in a pure state |9(] . It has been observed that the entanglement 
entropy of a block of L spins with the rest of the system at criticality follows a logarithmic 
scaling law in one-dimensional spin systems. The basic ingredient behind this result is that 
at criticality the correlation length diverges and the system becomes conformal invariant. 
This conformal symmetry leads to the logarithmic scaling law of the entanglement entropy 
at criticality 10|. When the criticality is induced by a quench the scaling behavior changes. 
The formation of defects as well as the entanglement entropy of a block of L spins with the 
rest of the system in transverse Ising model, when QPT is induced by a quench, has been 
studied in details by Cincio et al 

The behavior of the geometric phase in XY spin chain during a linear quench caused by a 
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gradual turning off of the magnetic field has been studied earlier In a recent note [12 
we have formulated a dynamical mechanism for QPT induced by a quench based on the 
fact that spin fluctuation in the vicinity of the critical point leads to the fluctuation of the 
Berry phase acquired by the ground state of the system when it evolves in a closed path. 
The probability of the formation of defects is determined by the two-point correlation of 
this geometric phase factor. We have studied the formation of defects in transverse Ising 
model during critical slowing down based on this formalism and the results are found to be 
identical with those derived by other authors js-sl. We have also studied the scaling law of 
the entanglement entropy in the transverse Ising model at criticality induced by a quench 
and it is shown that the prefactor now depends on the quench time and there is a restriction 



on the block size L 



13j. In this paper, we study the formation of defects as well as the 



scaling law of the entanglement entropy at criticality induced by a quench in some other 
one dimensional spin systems on the basis of this dynamical mechanism of QPT induced by 
a quench. 

The organization of the paper is as follows. In secHIlwe briefly review of our work related 
to the isotropic XY model (transverse Ising model) for comprehension of the formulation we 
have adopted. In sec lllll and sec JI VI we study the XX and XXX model respectively. In seclVl 
we deal with the Lipkin-Meshkov-Glick (LMG) model. We make a comparative analysis 
of the results obtained for different spin systems in sec JVIi Finally, we discuss results in 
seclVin 



II. THE XY MODEL 



The Hamiltonian for the XY model is given by 

^^^ = -11 (^^^'^i + ^^^r<i + ) • (1) 

Here A is the external magnetic field, 7 is the anisotropy parameter and erf (a = x,?/, z) are 
the Pauli matrices. The XY model with (7 7^ 0) falls into the free fermion universality class 
and is critical for A = 1. The isotropic XY model (7 = 1) corresponds to the transverse 
Ising model. We consider the Hamiltonian for the transverse Ising model 

^ = - E (^^'^^+1 + ^^^) ■ (2) 

i 
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For A >> 1 all the spins are aligned along the x-axis and the system is in the paramagnetic 
state. In the region A << 1 the system is in the ferromagnetic state with all spins either up 
or down. During critical slowing down we introduce a linear quench 

A(t<o) = --, (3) 

where Tg is the quench time. In the vicinity of the critical region, the nonadiabaticity induces 
spin fluctuation. This eventually causes fluctuation of the Berry phase acquired by a spin 
state when the ground state evolves in a closed path as the Berry phase acquired by a spin 



state is given by F = 
the quantization axis 
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cos6) where 6 is the angle of deviation of the spin axis from 



14j4l7l|. In the critical region, due to nonadiabaticity some spins get 
excited and defects (kinks) are formed. The Berry phase factor associated with such a spin 
is given by 

<Pko = ^ = 1(1 -cose,,), (4) 

where /cq is the momentum mode of the quasiparticle corresponding to the spin undergoing 
excitation and is the corresponding angle of deviation of the spin axis from the z-axis. 
In our formulation it is assumed that during criticality (pko undergoes stochastic fluctuation 
and it follows the simplest stochastic differential equation 

(i(/)fco(t) = -Uko(f>ko{'t)dt + drjit), (5) 

where is the frequency related to the energy e^g near criticality and r]{t) is a Gaussian 
white noise satisfying the moments 

{dvit)) = 0, 

{dr]{t)d7]{t')) =6{t- t')dt'. (6) 



This gives rise to the correlation 



(0fcoW) = O, 



(0.oW0;^o(O) = ^e--^o(*-*'). (7) 

At criticality during the quench time Tg, the spin states transit from the paramagnetic state 
to the ferromagnetic state so that the angle of deviation {Oko{Tg)) of a spin axis is it/2. The 
spin state undergoing excitation settles down at t = with the spin axis reversed so that 
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Okoi^) takes the value tt with certain probabihty pkfy This imphes that the random variable 
4>ko(t) acquires the value 1/2 and 1 at t = and t = respectively i.e, we have 4>ko{'Tq) = |; 
0fc(,(O) = 1. This suggests that the random variable 2(f)ko{Tq)(f)ko{0) attains the value 1 with 
probability pkg. Thus using eqn.(I7]) the excitation probability p^o can be written as 



The energy e^g near the critical point is given by 12 1 



k 

€ko = 2(1 - cos ko) = Asin'^Y ~ ^O' (9) 

for small ko- This gives 

p,^ = e-2-'=o%, (10) 
The number density of defects is given by 



'1 

This gives rise to the domain size corresponding to the Kibble-Zurek (KZ)correlation length 
^ ~ V^- This result is identical with that obtained by other authors 

At the critical point the ground states are characterized by the fact that the entanglement 
entropy 5*^ of a block of L spins with the rest of the system diverges like log L with a 
prefactor determined by the central charge c of the relevant conformal field theory [lo|. In 
fact, 

Sl^I I092L. (12) 



For the critical Ising model c = 1/2, 



Sl^]: log^L, (13) 




and the system falls into the free fermion universality class. It has been shown that the 
central charge c is related to the Berry phase factor 0, the phase being e*^'^'^ and just like 
the central charge c, the Berry phase factor also undergoes a renormalization group (RG) 
flow (is! such that < where L is a length scale. From this we have 

l^lx, = a In L = a log2L. (14) 

where a(< 0) is a parameter. 

For a pair of nearest neighbor spins (L = 2), a = |0| which corresponds to the entanglement 



of formation given by concurrence in a mixed state [?, 8|. So for the entropy of a block of L 
spins due to entanglement with the rest of the system in the pure state we can write 

Sl ~ 101 log2L. (15) 

At criticality the concurrence for the entanglement of a pair of nearest neighbor spins in 
transverse Ising model is given by \(f)\ = 0.18 [?! and so we find 

Sl ^ 0.18 log2L. (16) 

It is observed that the prefactor in eqn. fll3l) is very close to the value 0.18. We may mention 
here that the entanglement of a block of L spins with the rest of the system can be considered 
to be equivalent to the entanglement between a single spin representing the block spin with 
another spin represented by the rest of the system in block variable RG scheme. In view of 
this, Sl in eqn. flT^ can be considered as the concurrence C for the entanglement between 
the pair of this two block variable renormalized spins in a mixed state. The slight departure 
of the prefactor | in eqn.( ll3p from the value 0.18 in eqn.f ll6p may be associated with the 
block variable renormalization of the spin system which induces a change in the coupling 
constant. The factor = 0.926 is considered to be the correction factor associated with 

U. io 

the block spin variable. 

The introduction of a quench incorporates a new length scale given by the KZ correlation 
length ^ ~ So to evaluate the entanglement entropy immediately after the dynamical 
phase transition we have to rescale the Berry phase factor |0| in expression f|T5|) correspond- 
ing to the entanglement of a pair of nearest neighbor spins in terms of that of the block 
spins in the domain ^ with the rest of the system. This follows from the fact that after 
the formation of defects the system in the final state represents a kink-antikink chain with 



lattice space approximately given by ^. So we write 
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5(L,r,)=2Mi£^x 0.926 -3.7^^. (17) 
101 log2^ IriTg 

The factor 0.926 has been introduced as the correction factor of the block variable renor- 
malization as discussed above. Eqn.(17) shows that in the scaling law of the entanglement 
entropy there is a prefactor depending on This imposes a restriction on the maximum 
value of L, L^ax which is allowed in the system. As in the final state the system represents a 
kink-antikink chain with lattice constant ^, the maximum value of the entanglement entropy 



attained by the system corresponds to the entanglement of a block C, spins with the rest of 
the system along with the entropy associated with the entanglement of kink-antikink pair. 
It As we know for a bipartite spin 1/2 system the maximum value of entanglement entropy 
is 1, which at criticality is distributed over the whole chain, we can write for the entire 
chain jisl 

Smax = 2(0.18 log2^+ 1) X 0.926 = 0.12 In Tg + 1.85. (18) 
This suggests the relation ^^^'^''^ < 1 which implies IiqI 

In L< 0.03 {InTgf + 0.5 In Tg. (19) 

This restricts the block size L when QPT is induced by a quench and suggests that for 
sufficiently small Tg the block size L < ^. Thus we find that during critical slowing down 
due to the formation of defects, the scaling law for the entanglement entropy is valid only 
in a restrictive sense. 

We have compared our result with that obtained by Cincio et. al. jsj using standard tech- 



nique of density matrix analysis in [13|. Within the limited range of the block size L 
determined by eqn. (fT9l) the result was found to be in good agreement with that of Cincio 
et. al. 



III. THE XX MODEL 

The XX model given by the Hamiltonian 

^ = - E (^^'^i + + ^Y.^t (20) 

i i 

has two limit behavior. At A = 2 the system corresponds to a ferromagnetic state while at 
A = the system falls into the free boson universality class. The interval between these two 
points corresponds to the critical region. As the system transits through the point A = 2 in 
the vicinity of the critical point we introduce a quench so that the time dependent magnetic 
field is given by 

\{t < 0) = -2t/rg, (21) 

Tg being the quench time. At t = (A = 0) the system settles down in a free bosonic state 
when the ground state is ordered in the xy plane and we have pairs of spins with opposite 
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orientations. However, due to nonadiabaticity in the transition, this simple bosonic picture 
is disturbed and the system is characterized by domains having spins with same orientations 
i.e. defects (kinks) are formed. Following the analysis in the previous section, in the critical 
region we consider the fluctuation of the Berry phase acquired by a spin state when the 
ground state evolves in a closed path. Let the system initially be in the the region A > 2. 
As the system transits through the point A = 2, the ferromagnetic order is destroyed when 
certain spins are flipped. From eqn.(|l]) we can easily show that the Berry phase factor 
<t>ko{jq) takes the value 1 corresponding to the value Ok^i^Tq) = vr where ko is the momentum 
mode of the quasiparticle associated with the flipped spin. 

As time evolves, in the vicinity of A = due to nonadiabaticity certain spins get excited 
with probability p^j, and their orientations are reversed. Finally at A = 0, the system settles 
down with domains characterized by spins having this reversed orientation destroying the 
free bosonic property. This suggests that the random variable 0fco(O) attains the value 1 
corresponding to Ok^ = vr as follows from eqn.(j4l). From this we find that the random 
variable 4>ko{'Tq)4>ko{^) acquires the value 1 with probabihty p^o- Thus we have 



and from eqn.([7]) we find 



Pko = {<Pkoirq)(pkM). (22) 



Pko = ^e-'=o-^ (23) 



being the frequency associated with the energy in the critical region. 
Now we note that the energy of the quasiparticle associated with a spin state having the 
momentum mode ko is given by [2^ 

efc„ = A - 2 cos (24) 

being the number of sites. So taking the lattice constant a = ^ = 1, we find at A = 2, 
e^Q = 2(1 — cos ko) = 4 si'n?^ ~ k^ for small k^. This suggests that for we can write 

Pk, = ^e-'^^^'l (25) 

The number density of defects (kinks) formed when the system finally settles down at A = 
is given by 

n. = 5^ = l'- = (26) 



Incidentally it has the same scaling behavior with Tq as in the transverse Ising model and 
the KZ correlation legth is ~ 

To estimate the entanglement entropy at criticality we note that at A = the system 
belongs to the free boson universality class and the central charge of the relevant conformal 
field theory corresponds to c = 1. This gives the scaling law of the entanglement entropy at 
A = of a block of L spins with the rest of the system as 

~ i log^L. (27) 

As mentioned earlier the entanglement entropy of a block of L spins with the rest in a pure 
state is equivalent to the entanglement of formation of a pair of nearest neighbor spins in 
a mixed state given by concurrence C. As the concurrence C is associated with the Berry 
phase factor acquired by a spin state while evolving in a closed path we write C = 
where the phase is e^'^'"'^. Thus as in secUl] using a RG flow equation we can write for the 
entanglement entropy of a block of L spins with the rest of the system 

Sl ~ 101 log2L. (28) 

At the critical point A = 0, the system belongs to the free boson universality class, and the 
Berry phase factor |0| here is identical to the concurrence of a pair of nearest neighbor spins 
in an antiferromagnetic system which is given by C = |0| = 0.386 Now we take into 
account the correction factor 0.926 due to the block spin renormalization. In fact the Von 
Neumann entropy of a block of L spins with the rest of the system in pure state is here 
transcribed into the concurrence associated with the entanglement of a renormalized block 
spin with another spin corresponding to the rest of the system in a mixed state. This gives 
the effective value of |0| = 0.386 x 0.926 = 0.35 which is close to the value 1/3 derived from 
conformal field theory. Utilizing this result we now estimate the entanglement entropy of a 
block of L spins with the rest of the system at criticality near A = induced by a quench. 
Indeed in analogy to eqn. ffTTl) we write 



S(L, r„) = 4—^ — In L x 0.926 = 3.7^^. (29) 

101 In Tq In Tq 



This result is independent of the prefactor in eqn. (l28l) and thus is universal for all one 
dimensional spin systems. However, the maximum value of L, L^ax allowed in the system 
is different. Following eqn.( lT8|) . the maximum value of the entanglement entropy at A = 
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FIG. 1: In fig.l (left) the entanglement entropy of the XX model for various quench time (rg) is 



compared with the exact value plotted for A = from the result as reported in 
shows the variation of S{L,Tq)/ Smax with L/^^Tq. 



20]. Fig.l (right) 



is derived as 



Sn 



2(|(^| log2^ + l) X 0.926, 



0.25 In Tq + 1.85. 



(30) 



Now from the constraint S{L,Tq)/ Smax < 1 we obtain 



In L< 0.07 {In Tqf + 0.5 In Tq. 



(31) 



Thus we find that due to generation of defects, the scahng law for the entanglement entropy 
is restricted by the maximum value of the block size allowed in the system. This picture 
is in the vicinity of the point A = 0. However, in the region 2 > A > 0, the entropy of 
entanglement gradually decreases with the increase in the magnetic field following the same 
scaling law until at A = 2, the entanglement vanishes. 

It is pointed out by Latorre, Rico and Vidal 2l|] that the maximum entropy is reached 
at A = 0. The entropy is fitted by a logarithmic scaling law as given by eqnj ^7|l with 
a constant term which has been determined analytically by Jin and Korepin |22,]. These 
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results have been obtained without critical slowing down when nonadiabaticity comes into 
play. It is noted from eqn. fl29|) that when nonadiabaticity near criticality is introduced the 
entropy decreases with quench time. In figH] (left) we have compared our result with the 
exact value of the entropy at A = as reported in [2^ (where no quench was introduced) 
for Tq = 200, 400 and 800. As expected we have found that the entropy increases with 
decreasing quench time. At A = the entropy is maximum and we find that our result is in 
reasonable good agreement with the exact value for small value of Tq. Indeed, we find that 
for Tq = 200 it is close to the exact value and for higher values of Tq the entropy decreases. 
In figJU (right) we have plotted S{L,Tq) / Smax vs. Lj ^fTq. 



IV. THE XXX MODEL 



The Hamiltonian for the XXX model is given by 

^xxx = + + ^^^^1) + E (32) 

i i 

The critical behavior of this model is analogous to that of the XX model. It has two limit 
behavior. At A = 2 the system represents the ferromagnetic state and at A = the system 
corresponds to the antiferromagnetic state. The interval 2 > A > is gapless and hence 
critical. Let us denote the time evolution of the magnetic field by the relation 

A(t < 0) = -2t/Tq, (33) 

Tq being the quench time and analyze the dynamics of the system when it transits through 
the point A = 2. In the critical region, there is random fluctuation in the Berry phase and 
from the two point correlation of the Berry phase factor we derive the excitation probability. 
Actually, this is given by the same relation as in eqn. (l25l) and we write 

Pko = \e-^'o^^. (34) 

In this case the energy of the quasiparticle in the vicinity of the critical point is different 
from that of the XX model. 

A generalization of the XXX model is the anisotropic Heisenberg model which is known as 
the XXZ model given by the Hamiltonian 

i i 
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It may be mentioned that for A = the anisotropic system shows a gapless phase in the 
interval — 1 < A < 1. For A = 1 the system corresponds to the isotropic Heisenberg modeL 
For A > 2 the isotropic Heisenberg model corresponds to a ferromagnetic state. As the 
magnetic field decreases and passes through the critical point A = 2 some spins get flipped. 
For an excited state of the Heisenberg model belonging to the subspace of states where some 



spins are flipped the energy of the quasiparticle with momentum mode ko is given by 23| 



Cfco — 4(1 ~ cos ko) = 8sin^ ~ 2A;q, (36) 
for small ko- Thus from the relation (l34l) we can write 

and the number density of defects is given by 

^s = ^ r Pkodko = ^ r e~'-^^'hko = (38) 

_i 

Thus the number density of kinks scales like Tq ^ and the KZ correlation length is given by 

As the XXX chain at A = corresponds to the antiferromagnetic state, the behavior of the 
entanglement entropy for a block of L spins with the rest of the system is identical with that 
of the XX model. This is valid for S{L^ Tq) and Smax also. Just as in the XX model, in this 
model also, the entanglement entropy is maximum for A = which gradually diminishes as 
the magnetic field increases, and finally it vanishes at A = 2. 

V. THE LIPKIN-MESHKOV-GLICK MODEL 

The Lipkin-Meshkov-Glick (LMG) model demonstrates the mechanism of a phase transition 



for a many body system and was introduced almost forty five years ago [2J]. In contrast 
to the conventional spin models in the LMG model each spin interacts with all the spins 
of the system with the same coupling strength. This highly symmetric interaction pattern 
introduces the loss of the notion of geometry as there is no distance between the spins. 
This implies that we cannot consider the notion of a block of spins as a set of contiguous 
spins here. The symmetry of the Hamiltonian suggests that the ground state belongs to 
a symmetric subspace where all the spins are indistinguishable and this subspace restricts 
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the entanglement entropy of a block of L spins with the remaining spins. Interestingly 
the scaling behavior shows a similar pattern with that of the XX model where conformal 
symmetry plays a significant role at the critical region. But in the LMG model the scaling 
law seems to have nothing to do with any underlying conformal symmetry. 
In this note, we show that in case we introduce point-splitting regularization the regularized 
LMG Hamiltonian in the isotropic case (7 = 1) appears as a combination of Heisenberg 
(XXX) model and Ising model with magnetic field along the negative z-axis. In view of 
this, the regularized Hamiltonian follows the pattern of the usual spin systems where we can 
conceive of the notions of local interactions and translational invariance. The underlying 
conformal symmetry at the critical region of these systems is manifested in the LMG model 
and the similarity of the pattern of the scaling behavior of the entanglement entropy of 
this system can be understood. Thus the point-splitting regularization of the LMG model 
uncovers the hidden conformal symmetry at criticality which is lost when we take the sharp 
point limit. 

The Hamiltonian for the LMG model is given by 

^ = ^ E i^^^^ + + ^ E (39) 

i<j i 

N being the total number of spins. 

In what follows we shall consider the isotropic case with 7 = 1. The Hamiltonian (l39|l can 
be written in terms of the total spin operator 5*" = | ^ af as 

H = ^(^S^- {S'f - + 2A5^ (40) 

Let us consider that this total spin operator is located at a spatial point k. Now we introduce 
the point-splitting regularization so that we write 

S'^ = SkSk'Skk', (41) 

where k and k' are two adjacent sites with an infinitesimal distance k' — k = e. Thus S'^ 
corresponds to the product S^Sk' in the limit e — > 0. Taking 

Sk = lJ2^^^ ^." = ^E^"' (42) 
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where i and j are two adjacent sites with an infinitesimal distance. Considering only nearest 
neighbor interactions we can write the Hamiltonian fj40|) in the regularized form 



reg 



2 

N 



2Ar 

1 
2N 

-1. 



\ 1,j 1 / 



(43) 



In deriving eqn.f l43|) we have considered nearest neighbor interaction only because of the 
fact that the introduction of next nearest neighbor interaction gives rise to frustration in 
the antiferromagnetic phase of the XXX model which is unsolicited in the LMG model. In 
fact, this happens when the spin at the initial site interacts with all odd sites. To avoid this 
unwanted feature only the nearest neighbor interaction has been taken into account. 
Now considering 2^ as a dimensionless coupling constant J we write Hreg = Hi + H2 — I 
with 



Hi 



Ho 



J 



A V- 



■^(e-^m 



(44) 



(45) 



The Hamiltonian given by eqn. (143|) can be split into various other Hi and H2 such that both 
of them describe critical systems. The motivation behind this particular choice follows from 
the fact that cluster of even (odd) number of interacting spins represent bosonic (fermionic) 
systems. As the LMG model is characterized by the fact that each spin interacts with every 
other spin in the regularized Hamiltonian, we have chosen the critical systems such that one 
of these represents bosonic and the other fermionic features. In view of this we have taken 
Hi and H2 such that one represents the XXX model and the other Ising model. 
The criticality of this system can be studied by noting that this is a combination of the XXX 
model given by Hi and the Ising model with magnetic field along the negative z direction 
given by H2. It is observed that ^ denotes the intensity of the magnetic field. We know that 
the XXX model has two limit behavior. When ^ > 2, the XXX system is gapped and for 
A = 0, the magnetization is zero and the system is in an entangled (antiferromagnetic) state. 
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In the interval 2 > ^ > the system is gapless and this denotes the critical region. Noting 
that the dimensionless parameter J corresponds to being an integer we find that 

^ = 2 corresponds to the value |A| = j^. Since in an entangled state the minimum number 
of spins must be 2., i.e. Nmin = 2, at criticality |A| lies in the interval < |A| < 1. The value 
of A is independent of the real number of sites as we have treated the parameter 2^ as the 
coupling constant. Now it is observed that in the Ising chain given by H2 for |A| > 1 the 
system is in the ferromagnetic state when all spins are oriented along the negative z-axis. 
So in the interval < |A| < 1 as |A| is tuned from 1 to the spin system will undergo 
a transition when down spins will be excited and at A = all spins will settle down with 
opposite orientations. Evidently, with the evolution of |A| the spins cross a point when these 
are oriented along the x-axis. Hence, this corresponds to the double ferromagnet transition 
and in the interval < |A| < 1 spins evolve through a paramagnetic state. Thus during 
this transition in the interval < |A| < 1 the entanglement entropy of the spins evolves 
through a nonzero value which is similar to that of the transverse Ising model. In the 
regularized Hamiltonian we can consider the time dependent magnetic field corresponding 
to the summation of that given by eqns.(l3]) and ( |33l) and so we write 



A(t < 0) = + - = --, (46) 



Tq 'Tq Tq 



Tq being the quench time. We have taken ^ for the Ising system with positive sign as the 
magnetic field here is in the negative z-direction. Finally, when the system settles down at 
if: = (A = 0), the number density of defects is estimated and is given by = rii + 
where ni {113) is given by eqns.l lTT]) (( l38l) ) in the thermodynamic limit. From this we obtain 

11 11 2^2 + 1 
n4 = m + ng = = + — = — — — . (47) 

_1 

Thus we note that in LMG model also the number density of kinks formed scales as ^ . 
However, in this model a correlation length characterizing the typical distance between 
defects cannot be introduced though one can estimate the fraction of fiipped spins after the 
quench. 

The entanglement entropy in the LMG model of a block of L spins with respect to the rest 
of the spins can be derived in an analogous way. The entanglement entropy in the critical 
region < |A| < 1 is given by the summation of the entanglement entropy of the XXX 
model and the transverse Ising model at criticality. As the entanglement entropy for a block 
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of L spins with the rest of the system is given by Si{L) ~ | log2L for the XXX model and 
that for the transverse Ising model is given by 5*2 (L) ~ | log2L, the entanglement entropy 
in the critical region < |A| < 1 for the LMG model corresponds to the scaling law in the 



thermodynamical limit as 



S'(L) ~ ^ log2L + ^ log2L ~ ^ log2L, (48) 



which is identical with that obtained by Latorre, Orus, Rico and Vidal [25j. One should note 
that for A > 1, the system represents a polarized product state. Thus when the LMG model 
is recast in the regularized form the critical region manifests the same logarithmic scaling 
law as observed in XXX model and the transverse Ising model which are characterized by 
conformal symmetry at the critical region in these systems. Though the LMG model appears 
to have nothing to do with conformal symmetry at the critical region, the point-splitting 
regular izat ion reveals an underlying conformal symmetry at criticality in this system which 
is lost in the sharp point limit. 

Now we consider that the criticality is induced by a linear quench f H6|) . As the regularized 
Hamiltonian suggests that at the critical region this corresponds to a combination of the 



XXX model and the transverse Ising model so following eqns. f[T7j ) and fl29|) . we write 

S(L.r,) = i|i^*£x 0,926. (49) 

101 log2^ 

where |0| = + \4>\ising = 0.386 + 0.18 = 0.566 and ^ (~ is the associated KZ 

correlation length. The final result is independent of |0| and we obtain the universal value 

In T 

S{L,rq) = 3.7^. (50) 

In Tg 

Here the maximum value of L, Lmax is different from that of the XXX and the transverse 
Ising model. In analogy to eqns. (fTSj l and (l30l) we find that the maximum value of the 
entanglement entropy 

SmaxliMG = 2(101 log2^+ 1) X 0.926 ^ 0.36 In t, + 1.85. (51) 

So from the constraint S{L,Tq)/ Smax < 1 we find 

In L\lmg < 0.097 {In r^f + 0.5 In r^. (52) 
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In the sharp point hmit the LMG model does not allow any correlation length. However, 
we can consider the coherence number 26| which can be identified with the KZ correlation 
length ^ in the regularized formalism. It is noted that in the critical region < A < 1 the 
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FIG. 2: In fig. 2 (left) we have compared our result for S{L,Tq) for Tq = 200, 400 and 800 with the 
values derived from that reported in [20] at A = extrapolating it to the thermodynamic limit. In 
fig. 2 (right) we have plotted our results for L = 16, 32 and 64 in the region Tg > 200. 



entanglement entropy is maximum at A = and decreases with the increase in the magnetic 
field until at A = 1 it vanishes when the system represents a product state. However, within 
this region the scaling behavior remains the same. 

It may be pointed out here that recently the LMG model has been studied in which the 

system is dragged adiabatically through the critical point 2^]. Using the Landau-Zener 

_ 3 

formula [28] the fraction of the fiipped spins has been found to scale like ^ , being the 

_ 1 

quench time. The failure of obtaining the scaling behavior Tq ^ as derived here may be 
related to the definition of the defect density in this formalism. Indeed, here the degree of 
adiabaticity is estimated through the residual energy Eres given by Eres = -E'/m — Egs where 
Efin is the average energy in the final time-evolved state and Egg is the ground state energy. 

It is found that the system has three regions: (i) for fast quenches Eres persists independent 

_ 3 

of Tq, (ii) for slower quenches Eres persists ~ Tq ^ and (iii) for further slowing Eres persists 



~ T,. 



It is noted that for fast quench, the transition is nonadiabatic. However, for slow 



transitions the adiabatic dynamics is quite different from our present formalism. In fact, 
the result we have obtained here is fully determined by the nonadiabaticity condition in the 
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FIG. 3: The entanglement entropy S{L,Tg) vs. Lj ^Pf^. 



vicinity of the critical point. 
In fig|2] (left) we have compared our resu 
with the value for A = derived in 



t for the entanglement entropy for the LMG model 
20| extrapolating it to the thermodynamic limit. It 
is noted that when we introduce a quench the entropy decreases with the increase in the 
quench time r^. So it is expected that the value reported in [201] will be in agreement with 
our result for small r^. However for very small r^, the maximum value of the block size 
L will be very small which follows from eqn. (l52l) . In view of this we have taken moderate 



values of r„ = 200, 400 and 800 for comparison. It is noted that for r„ = 200 the result is 



very close to the value derived from 20 1. 

Caneva tt. al. [!27i] have computed the entanglement entropy at A = for various values 
of Tg when the system is adiabatically drugged through the critical point for finite system 
size. According to their result, for fast quench — )■ the entanglement entropy tends to 
be vanishing and for slow dynamics Tg ^ oo the entropy picks up the value it assumes 
in the final ground state Sgs (A = 0) = 1 independent of the system size. Between this 
two limiting behaviors the entropy reaches a size dependent maximum at Tg nearing 10 and 
then decreases with the increase in Tg. In our formalism we find that for Tg nearing 10, the 
maximum value of the block size L will be very small as follows from eqn.f l52p . In view of 
this we have taken moderately large Tg > 200 for comparison with their result. Caneva et. 
al 271] have computed their result for block size L = N/2, N being the system size. However 
our result has been computed in the thermodynamic limit. In figj2] (right) we have plotted 
our results for L = 16, 32 and 64 and Tg > 200. We find that our result is in reasonable good 
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agreement with that of Caneva et. al. [27|. In fig|3]we have plotted S{L,Tq) vs. L/^/t^. 



VI. DEFECT FORMATION AND MAXIMAL ENTANGLEMENT ENTROPY IN 
DIFFERENT ONE DIMENSIONAL SYSTEMS: A COMPARATIVE STUDY 

We have studied here defect formation and the scahng of entanglement entropy in QPT 
in different one dimensional systems induced by a quench. It is found that they all show 



5 
4.5 
4 




1.5 

n 

3 

2.5 

2 








3 4 5 6 7 8 
Inx 

q 



FIG. 4: Fig. 4 (left) shows the variation of the number density of defects ni, n2, n^, n4, with 
quench time Tq in isotropic XY, XX, XXX and LMG models respectively for different values of Tq. 
Fig. 4 (right) shows the maximum value of the entanglement entropy Smax for different values of 
Tq in isotropic XY [Smaxx)-, XX (XXX) {Smax2) ^iid LMG {Srnaxs) models. It is noted that Smax 
is identical for XX and XXX systems. 



a universal behavior. The number density of defects in all these systems scales like Tq ^ 
with difference only in the prefactors. Also the entanglement entropy at criticality takes 
a universal value having the same scaling pattern with block spin size L. However, there 
is a restriction on the maximum value of L allowed in different systems. The maximum 
value of the entropy in all these systems follows the same pattern with difference only in 
the coefficients. In all these systems the constraint on the maximum value of the block size 
L restricts the validity of the scaling law for any arbitrary L. In fact, the main implication 
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of the scaling law for the entanglement entropy is that with the addition of every spin, 
entanglement increases. However, in QPT induced by a quench this is not valid beyond 
the upper limit of the value of L allowed in various systems. Indeed, the nonadiabatic 
transition due to quench makes the scaling law a restrictive one and cannot be treated in 
the conventional sense. 

In figJH (left) we have plotted the number density of defects formed during QPT induced 
by a quench in different spin systems for different values of Tq. The density of defects for 
the isotropic XY, XX, XXX, LMG models are denoted by ni, n2, na, respectively. In 
figJH (right) we have plotted the maximum value of entanglement entropy S„iax for different 
values of Tq in isotropic XY, XX and LMG models. Smax is identical for XX and XXX 
systems. 



VII. DISCUSSION 



We have analyzed here the defect formation and the scaling law of the entanglement 
entropy in QPT in several one- dimensional systems at critical slowing down when nonadi- 
abaticity plays a dominant role. We have argued that in the vicinity of the critical point 
nonadiabaticity in QPT causes spin fluctuation which in turn makes the Berry phase factor 
a random one. From the two-point correlation of the Berry phase factor we have estimated 
the probability of the generation of defects. Indeed, in a recent paper [l2| we have esti- 
mated the number density of defects as well as the spin-spin correlation at criticality of 
the transverse Ising model from this dynamical mechanism of QPT induced by a quench. 
The results are found to be identical with those derived from the standard Landau-Zener 
transition probability studied by other authors js-Sj. Here we have analyzed several other 
one dimensional spin systems using this dynamical mechanism. Our analysis suggests that 
this formalism represents the universal dynamics of QPT induced by a quench. 
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